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Abstract 

We study the neutralino relic density and the neutralino-proton elastic 
scattering cross sections in the minimal supersymmetric standard model 
(MSSM) with CP-violating phases. We include all the final states to the 
neutralino pair annihilation cross section at the tree level, taking into ac- 
count the mixing between the CP-even and CP-odd neutral Higgs fields. 
We show that variations of the relic density and the elastic scattering cross 
sections with the CP- violating phases are significant. In particular, inter- 
ferences between the major contributions in the neutralino annihilation 
cross section can cause considerable enhancement of the relic density. 
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1 Introduction 



Recent progress in cosmological observations has much impact on determination of 
various cosmological parameters. The analyses of the anisotropy in the cosmic mi- 
crowave background radiation using the first year WMAP data provide the relic abun- 
dance of a cold dark matter (CDM) as [1,2] 

Qcdm/i' = 0.1126l°:°°«, (1) 

where i^cDM is the CDM energy density normalized by the critical density and h ~ 
0.7 is a parameter in the Hubble constant Hq = 100 h km/sec/Mpc [3]. The current 
allowed region of the CDM relic density is expected to be narrowed by the analyses 
of increasing WMAP data and forthcoming data from the future project Planck [4]. 

One of the most promising candidates for the CDM is the lightest supcrparticle 
(LSP) in the minimal supcrsymmctric standard model (MSSM) with R-parity conser- 
vation [5]. The R-parity guarantees that the LSP is stable. In this model, the LSP is 
typically the lightest neutralino [6] which is a linear combination of neutral gauginos 
and higgsinos 

X = N^,B + N^2W^ + N^,H', + NuHl (2) 

where B is the f/(l)y gaugino (bino), and is the neutral SU{2)l gaugino (wino). 
ifj* and iJ^ the two neutral higgsinos with opposite hypercharges. The coefficients 
Nii {i = 1,2,3,4) are the elements of the 4x4 unitary matrix N which diagonalizes 
the neutralino mass matrix [5]. 

The lightest neutralino has weak but finite couplings to the quarks. This implies 
that it may be possible to detect the neutralino CDM directly in a detector. There 
have been various ongoing and planned activities on the direct detection of Weakly 
Interacting Massive Particles (WIMPs) [7-10]. Current experimental upper bounds on 
the elastic scattering cross section of the WIMP off proton are around 10^® pb and 1 
pb for spin-independent (SI) interactions and spin-dependent (SD) ones, respectively. 
Future experiments are expected to drastically improve the sensitivity. Thus we may 
observe a first signal of a supcrparticle in the WIMP searches. 

There have been a lot of theoretical analyses on the direct detection [11-16] and 
the relic density [17-27] of the neutralino CDM in the MSSM [28]. However, be- 
cause of a large number of free parameters in the MSSM, most analyses assume some 



2 



pattern of the supersymmetric parameters to reduce a number of free parameters, 
such as a common sfermion mass, a definite relation on the gaugino masses, ignoring 
generation mixing or complex phases in the supersymmetric parameters. Considering 
increasing sensitivities of experimental projects, it is now important to relax these 
conditions to examine the possible variations of the CDM observables under more 
general conditions. 

In particular, the dependence of the CDM observables on the CP- violating complex 
phases needs to be carefully examined. It has been realized that non-zero phases of 
the trilinear scalar couplings associated with the third generation sfermions {At, A^, 
At-) or the Higgs mixing mass parameters (fi) can induce the mixing between the CP- 
even and CP-odd neutral Higgs fields [29,30]. In Ref. [31], the direct detection cross 
sections of the lightest neutralino in the presence of CP violation were studied ignoring 
this mixing. An analysis of the direct detection cross sections with supersymmetric 
CP violation taking into account the CP even-odd mixing in the neutral Higgs sector 
has been done in Ref. [32]. It was shown that the CP- violating phases can reduce or 
enhance the neutralino-nucleus cross sections significantly. 

On the other hand, a full analysis of the neutralino relic density Q^h"^ including 
all the contributions to the annihilation cross section in the MSSM with CP violation 
is missing. The effect of CP violation on the neutralino annihilation cross section into 
the fermion pairs through the sfermion exchange was studied in Ref. [33]. It was shown 
that a non-zero phase in the squark left-right mixing removes the p-wave suppression 
of the annihilation cross section for the fermion pair final states, and greatly enhances 
the cross section. In Ref. [34] , the effects of a CP-violaing phase in the trilinear scalar 
couplings on the neutralino pair annihilation cross section were examined, including 
all the tree level contributions and taking into account the mixing between the CP- 
even and CP-odd neutral Higgs fields. It was found that, the annihilation cross section 
can be enhanced typically by factors of one to four; the enhancement can be huge near 
the resonances where the relic density increases by a large resonant contribution of the 
s-channel Higgs exchange. Recently, the dependence of Q^h'^ on CP-violating phases 
through supersymmetric loop corrections to the bottom-quark mass was examined, 
taking into account the mixing between the CP-even and CP-odd neutral Higgs fields 
and focussing on the sfermion exchange and the Higgs boson exchange contributions 
to the fermion pair production [35]. It was shown that, near the resonances, the relic 
density is sensitive to CP-violating phases mainly through the loop corrections to the 
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bottom-quark mass. 

In the present work, we re-investigate the effects of CP-violating complex phases of 
supersymmetric parameters on the neutralino relic density and the neutralino-proton 
elastic scattering cross sections. We include all the contributions to the neutralino 
annihilation cross section at the tree level, taking into account the mixing between 
the CP-even and CP-odd Higgs fields. We examine the effects of CP violation outside 
the resonance regions. 

This paper is organized as follows. In section 2, we declare our assumption on the 
parameters in the MSSM with supersymmetric CP violation, and briefly describe the 
structure of the relevant interactions. In section 3, we give the formulae on the direct 
detection cross sections. In section 4, we briefly sketch our method to compute the 
relic density of the lightest neutralino. In section 5, we present our numerical results. 
Finally we give conclusions in Section 6. 

2 The MSSM with supersymmetric CP violation 

In general, the MSSM has a lot of free parameters [5]. In the prerent work, we examine 
the effects of the following parameters 

Ml, M2, M3, mp Af, tan/3, m^, (3) 

where Mi, M2 and M3 are the mass parameters for the bino, the wino and the gluino, 
respectively. // is the Higgs mixing mass parameter, is the supersymmetry break- 
ing common mass parameter for the sfermions. Af is the common trilinear scalar 
coupling for the third generation, while the ones for the first two generations are ne- 
glected, tan P is the ratio of the vacuum expectation values of the two neutral Higgs 
fields, rrij^ is a parameter which coincides with the CP-odd Higgs mass in the absence 
of CP violation [29,30]. In eq.(3), Mj (i — 1,2,3), /j, and Af can have CP-violating 
phases in general. 

For the gaugino masses, we consider two cases: grand unified theory-like (GUT- 
like) case and anomaly-mediated supersymmetry breaking-like (AMSB-like) case. For 
the GUT-like case, the absolute values of the gaugino masses are related as 

|Mi| = Itan^^p^lMsl, IM3I = — IM2I, (4) 
3 a2 
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where 6w is the Weinberg angle, and ^2 and are the gauge couphng constants for 
SU{2)l and SU{3)c gauge groups, respectively. In the GUT-like case, the neutralino 
LSP can be bino-like (|Mi| < or higgsino-like (|Mi| > For the AMSB-like 
case, they are related as [36] 

\Mi\ = 2.8IM2I, IM3I = -8.3IM2I. (5) 

In this case, the neutrahno LSP can be wino-hke (IM2I < or higgsino-like (IM2I > 
Note that the SU(5) grand unification or the anomaly- mediated supersymmetry 
breaking scenario implies the similar relations not for the absolute values |Mj| {i — 
1, 2, 3) but for Mj themselves. 

We take a convention that M2 is real. Also, Af and M3 are assumed to be real. 
The CP- violating phases we examine in the present work are the phases of Mi and /i 

Ml = |Mi| exp(i6'Mi), A* = exp(i6'^). (6) 

It has been realized that CP- violating phases in Af or /i can induce the mixing 
between the CP-even and CP-odd Higgs fields through radiative corrections [29,30]. 
The MSSM contains two CP-even neutral Higgses 0i and 02 with hypercharges Y{(f)i) 
— — y(02) = ~l/2 and a CP-odd neutral Higgs a as physical fields. In general, non- 
zero phases in Af or /i induce (pi-a and 02-a mixings at the one-loop level, so that the 
Higgs boson mass-squared matrix M.^ becomes 3x3. We calculate the Higgs mass- 
squared matrix using the one-loop effective potential, including the contributions of 
the third generation fermions and sfermions. The mass eigenstates (r = 1, 2, 3) 
are related with the CP eigenstates (0i,02,a) by a 3x3 rotation matrix Oh as follows 




The Higgs boson mass eigenvalues are obtained as OhM.\OJj = diag(m^(,,m^(,,m^(,). 
Note that the parameter in eq.(3) is not a mass eigenvalue in general. The CP 
even-odd mixing in the neutral Higgs sector implies (0//)j3 7^ (i = 1, 2) so that 
interactions of the neutral Higgs bosons are significantly modified [29,30]. 

If the mass parameters Mi and fi have non-zero phases, the elements of the neu- 
tralino mixing matrix N have an imaginary part in general. In this case, the interac- 
tions of the neutral Higgses with the lightest neutralino and the fermions (/) relevant 
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for the present analysis have the following structure 

^xxHO = j:x{Cr''°-Cf^'''l.)xH'^, (8) 

r=l 

^ffHO = E tfiC'/'''-C'/"'^.)fH'r- (9) 
f=u,d,--- r=l 

Without CP violation, either the scalar coupling Cs or the pseudoscalar coupling 
Cp is vanishing for every mass eigenstate [5]. In the presence of CP violation, 
however, both Cs and Cp are non-zero for every mass eigenstate. The neutralino- 
fermion-sfermion interactions have similar structure 

%x -EE xiC^s"' - cf^'l^)f~U + h.c. . (10) 

/=u,d,---i=l,2 

We have neglected generation mixing. The two mass eigenstatcs fi {i = 1,2) for the 
sfermions are given by linear combinations of the gauge eigenstatcs fi and fji with 
the same flavor. The expressions for the coupling constants in equations (8), (9) and 
(10) can be found in Ref [24,32]. 

3 Direct detection of the neutralino dark matter 

In this section, wc give the formulae for the clastic scattering cross section of the 
lightest neutralino off proton XP XP- The fundamental process for this scattering 
is the neutralino-quark scattering XQ XQ- the non-relativistic limit, the relevant 
effective four-fermion interaction Lagrangian reads [13] 

^xq = dq{xYl5X){Ql,^l5q) + fq{Xx){Q(l)- (11) 

The couphng for the spin-dependent (SD) interaction consists of the t-channel Z boson 
exchange contribution and the s-channel squark exchange contribution [13] 

where T^g represents the isospin charge of the quark, nig, m^., and niw denote 
the masses for the quark, the squark, the LSP and the W boson, respectively. The 
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coupling for the spin-independent (SI) interaction contains the t-channel neutral Higgs 
boson exchange contribution and the s-channel squark exchange contribution 

3 r'XxH^f^qqH'^ . \r<xmi\2 \r'xmi\2 

JQ — 2^ ^2 A ^2 \2- y^'^) 



r=l "^i/o j=i 2 '"'Qi 

In general, the effective Lagrangian of neutralino-quark interaction (11) contains other 

terms: (xO'''75X)(^7m5)> (xx)(^75g), (X75X)(w) and (x75X)(^75g)- However, in the 
non-relativistic limit, these operators are negligible compared to the contributions of 
eq. (11), since contributions of these operators are suppressed by the small velocity 
of the neutralino and the proton [12,31]. We ignore them in the present analysis. 

One has to convert the quark level Lagrangian (11) to the hadronic effective La- 
grangian. The neutralino-proton SD interaction can be found to be described by the 
following effective Lagrangian: 

= dp{xri,x)iPWP). (14) 
The expression for the coupling constant dp is 

dp ^ i: d,A^\ (15) 

q=u,d,s 

Here the parameters A^^ are defined by (p|?7/^75g|p) = 2snA'^\ where is the spin 
of the proton. In the present analysis, we take [37] 

A(f) = 0.74, aJ'^ = -0.4, A^f = -0.12. (16) 

The SD cross section for the elastic scattering of the neutralino can be written as 



a 



SD _ ^_2j2 
XV ^ 



(17) 



where = m^rrip/ {m^ + nip) is the reduced mass. 

The neutrahno-proton SI interaction can be described by the following effective 
Lagrangian: 



4p = fpixx)iPp)- (18) 



The coupling constant fp is given by [38, 39] 



f q=u,a,s y q=c,o,t 'J 
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where f^Q = 1 — J2q=u,d,s frq ^ ^^^d the parameters f^^ are defined by ip\mqqq\p) ~ 
nipf^q {q = u, d, s). In our numerical analysis, we take [40] 

= 0.020, /S = 0.026, /i^i^ 0.118. (20) 

We do not include the contributions of twist-2 operators [13], since they are typically 
sub dominant. The SI cross section for the elastic scattering of the neutralino can be 
expressed as 

= ^"'m- (21) 

4 Relic density of the neutralino 

In this section, we briefly describe our method to compute the neutralino relic density 
in the MSSM with supersymmetric CP violation [28,41]. 

We need to evaluate the relic density at present, starting from thermal equilibrium 
in the early universe. The time evolution of the neutralino number density in the 
expanding universe is described by the Boltzmann equation 

^ + ^Hn^ = -KxWK-W']' (22) 

where H is the Hubble expansion rate, and is the number density which the 
neutralino would have in thermal equilibrium. Information of the MSSM Lagrangian 
is encoded in the cross section o"^^ of the neutralino pair annihilation into the standard 
model particles. The quantity {crxx''^M0i) represents the thermal average of crj^;^i'M0b 
where vm01 is a so-called M0ller velocity which can be identifled with the relative 
velocity between the two colliding neutralinos. 

In the early universe, the neutralino is assumed to be in thermal equilibrium where 
— n^. As the universe expands, the neutralino annihilation process freezes out, 
and after that the number of the neutralinos in a comoving volume remains constant. 
Using an approximate solution to the Boltzmann equation (22), the relic density 
— m^n^ at present is given by 

1.66 (TA^ 3 1 

where x = T/m^ is a temperature of the neutralino normalized by its mass. 
and are the present temperatures of the neutralino and the photon, respectively. 
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The suppression factor (T^/T^)^ ^ 1/20 follows from the entropy conservation in a 
comoving volume [42]. Mpi denotes the Planck mass, xp is the value of a; at freeze-out, 
and is obtained by solving the following equation iteratively: 



where is the Newton's constant, and (~ 81) represents the effective number of 
degrees of freedom at freeze-out. 

Accurate calculation of the relic density requires a careful treatment of the thermal 
average in eq. (22) and the annihilation cross section in it. In literatures, expansion of 
the thermal average in powers of the temperature {(t^^vm0i) ~ a + 6x is widely used. 
However, it is known that the expansion method causes large error when cr^^ varies 
rapidly with the energy of the neutralinos, hence, in gerenal, one has to use the exact 
thermal average written as an integration with respect to a Mandelstam variable [21] 



where Ki {i — 1,2) are the modified Bessel functions. The cross section cr-^xi^) ^ 
complicated function of s in general, so we have to evaluate the above integration 
numerically to calculate the thermal average. 

In the annihilation cross section a^x in eq. (22), there are a lot of final states: 



XX fl {r.s = 1,2,3), W+W' , ZZ, W^H^, ZH^ (r = 1,2,3). 



contributions. However, depending on the parameters, the other final states can also 
play a considerable role [20,22,24]. 

The expressions for the full cross section without CP violation at the tree level are 
found in Ref. [26]. We extend the analysis of Ref. [26] to incorporate the effects of CP 
violation. We have derived the full expressions for the neutralino annihilation cross 
section with CP- violating phases at the tree level, taking into account the modified 
interactions (8) and (9) in the presence of the mixing of the CP-even and CP-odd 
neutral Higgs fields. In the present work, we shall present our numerical results in 
the following. The analytic expressions for the annihilation cross section and more 
comprehensive analyses of the relic density in the presence of CP violation will be 
presented elsewhere. 




(24) 




(25) 



Among these final states, the fermion pair final states // often give the dominant 
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In the present analyses, we neglect coannihilation effects [43-49], though they are 
crucial when the LSP is higgsino-like or wino-lilce. The investigation of the effect of 
CP violation on the coannihilation cross sections is left for future work. 

5 Numerical results 

In this section, we present our numerical results. We examine the sensitivity of the 
SI and SD cross sections and the relic density on the phases in eq.(6). We search for 
the parameter sets allowed by the la range of the WMAP constraint 

0.0946 < Q^/i^ < 0.1286, (26) 

and present typical results in four cases: bino-like LSP case, mixed LSP (bino-higgsino 
mixing) case, higgsino-like LSP case and wino-like LSP case [50]. 

Let us start with showing the effects of the phase Qu^ . One of the typical results for 
a bino-like LSP case is illustrated in Figures l(a)-(d). In these figures, we take = 
300 GeV, Af = 300 GeV, /i = 200 GeV and M2 = 150 GeV with the GUT-like gaugino 
mass relations (4). The parameter rrij^ is fixed at rrij^ = 500 GeV for all the figures 
in the present work. The solid line, the dashed line and the dotted line correspond 
to tan/9 = 10, 30 and 50, respectively. In Fig. 1(a), we plot the SI neutralino-proton 
cross section cr^p in eq.(21) as a function of 6mi- There occur cancellations in the SI 
cross section for non- vanishing values of [31,32]. 

Fig. 1(b) shows the dependence of SD cross section in eq.(17) on 6mi- The 
dependence is milder than that in the SI case, but it shows a sizable variation of about 
a factor three. The case = tt gives the smallest value of a^^. 

In Fig. 1(c), we plot the behavior of the neutralino relic density Qy.h'^. The region 
between the two dash-dotted lines in this figure corresponds to the 2a range of the 
WMAP constraint (26). For tan/3 = 30 and 50, the relic density shows significant 
dependence on the phase 6mi- In particular, we see that, for tan/3 = 50, CP violation 
is necessary to find the relic density allowed by the WMAP constraint. The relic 
density is reduced for larger tan f3, since the annihilation cross section is enhanced via 
larger coupling of the neutralino with bottom-quark and bottom-squark if the LSP 
includes non-zero higgsino components. 

The relic density is enhanced around 6mi ~ 27r/3, 47r/3 for tan/3 = 30 and 50 
in Fig. l-(c). This enhancement results from the interferences between the major 
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contributions in the ncutralino annihilation cross section. For Om^ = 0, the dominant 
contribution to the annihilation cross section a^^, comes from the bottom-quark pair 
creation XX ~^ bb. In this process, a bottom-squark exchange diagram and a Higgs 
exchange diagram contribute constructively for 6mi = 0, and this leads to a too small 
relic density for tan/3 = 50. On the other hand, for 6mi ~ 47r/3 in the same 

figure, these diagrams contribute destructively. This interference leads to suppression 
of the annihilation cross section, and the relic density is enhanced to fit into the 
WMAP-allowed region for tan f3 = 50. 

The phase is constrained by experimentel constraints on the EDMs of the 
electron, the neutron [51-53] and ^^^Hg [54]. In the present work, we do not perform 
a complete analysis of the EDM constraints, but only present results for the electron 
EDM. In Fig. 1(d), we show the absolute value of the electron EDM de- The current 
experimental bound for the electron EDM is \de\ < 1.6 x lO^^^e ■ cm [55]. The upper 
bound is shown by the dash-dotted line in the same figure. With this parameter set, 
\de\ is beyond the upper bound for 6mi 7^ 0, vr. To satisfy the EDM constraint for 6mi 
7^ 0, TT, we have to take heavier mass parameters. This typically results in a too large 
relic density compared to the WMAP constraint. 

Figures 2(a)-(d) correspond to a mixed LSP (bino-higgsino mixing) case. In these 
figures, we take niA = 500 GeV, rrif = 800 GeV, Aj = 800 GeV, = 200 GeV and Ma 
= 300 GeV with the GUT-like gaugino mass relations (4). The SI cross section shows 
cancellations at non-zero values of for tan/3 = 30 and 50. The SD cross section is 
larger than the bino-like case, since for the mixed LSP case with bino-higgsino mixing, 
either A^i3 or N14 is sizable so that the Z boson exchange contribution in eq.(12) is 
increased. The SD cross section shows mild dependence on 9mi- The relic density in 
Fig. 2(c) again depends on 6mi significantly, and WMAP-allowed regions appear for 
non-zero 9mi in the case of tan (3 = 10 and 30. The WMAP-allowed regions for tan/3 
= 10 also satisfy the EDM constraint as seen in Fig. 2(d). 

For a higgsino-like LSP, the cross section a^^ is relatively large, and the relic den- 
sity is typically too small to satisfy the WMAP constraint. In order to find WMAP- 
allowed regions, we have to choose relatively large mass parameters to enhance the 
relic density. Figures 3(a)-(d) correspond to a higgsino-like LSP case. In these fig- 
ures, we take = 500 GeV, m^- = 1 TeV, Af = 500 GeV, /i = 700 GeV and M2 = 
2.5 TeV with the GUT-like gaugino mass relations (4). In this case, cancellations in 
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the SI cross section 0"^^ do not occur, but 0"^^ still depends on Om^ as seen in Fig. 3(a). 
The SD cross section in Fig. 3(b) is suppressed because |A^i3| is nearly equal to |A''i4| 
for m\A/ ^ < |Mi|, M2 so that the Z boson exchange contribution in eq.(12) is 
small. The relic density in Fig. 3(c) is not sensitive to 9mi- In the pure higgsino 
limit, the mass and couplings of the LSP do not include 6'mi at all. This explains the 
insensitivity to Omi- 

Figures 4(a)- (d) correspond to a wino-like LSP case. In these figures, we take 
mA = 500 GeV, = 2TeV, Af = 500 GeV, /x = 2TeV and M2 = 1.8 TeV with 
the AMSB-like gaugino mass relations (5). For a wino-like LSP, the relic density 
is typically smaller than the one in a higgsino-like case so that we have to choose 
relatively large mass parameters to find WMAP-allowed regions. We see that the 
cross sections a^p, and are not sensitive to 9mi in this case. 

The 9^ dependence of cr^p, a^, Q^h"^ and \de\ is shown in Figures 5(a)-(d), 6(a)- 
(d), 7(a)-(d) and 8(a)-(d) for the same choice of parameters as in Figures l(a)-(d), 
2(a)-(d), 3(a)-(d) and 4(a)-(d), respectively. We see that the effects of 6*^ are also 
relevant. Cancellations can occur in the SI cross section for non-trivial values of 9mi- 
The SD cross section is less sensitive to 9^ compared to 9mi- Variations of the relic 
density are again significant for the bino-like LSP and the mixed LSP, while the relic 
density is not sensitive to for the higgsino-like LSP and the wino-like LSP. Note that 
the EDM constraint is much severer for 9^ than 9mi, so that a sizable CP- violating 
phase 9^ (7^ 0, vr) is experimentally ruled out. 

In all the figures, the parameter sets do not correspond to resonance regions where 
the relic density is significantly reduced by resonant annihilation through the s-channel 
Higgs poles. By choosing an appropriate value of m^, we can hit the resonance regions 
where Q^h"^ is much smaller. Note that, near the resonances, the dependence of Q^h'^ 
on CP-violating phases through supersymmetric loop corrections to the bottom-quark 
mass is significant for large tan f3 [35] . 

In the present analysis, we have neglected coannihilation effects [43-49]. For Fig. 
1(c) and Fig. 2(c), the coannihilation effects are really negligible. However, for 
Fig. 3(c) and Fig. 4(c), there exist a lighter chargino and/or a next-to-lightest 
neutralino which are nearly degenerated with the LSP. This implies that inclusion of 
coannihilation effects is essential for Fig. 3(c) and Fig. 4(c), and this will result in a 
too small relic density compared to the WMAP constraint (1). The investigation of 
the effect of CP violation on the coannihilation cross sections is left for future work. 
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6 Conclusions 

We have studied the neutrahno rehc density and the neutrahno-proton elastic scatter- 
ing cross sections in the MSSM with CP-violating phases. We have included all the 
final states to the neutralino pair annihilation cross section at the tree level, taking 
into account the mixing between the CP-even and CP-odd Higgs fields. We have 
demonstrated that the variations of the relic density and the elastic scattering cross 
sections with the CP- violating phases are significant. 
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Figure 1: The results for = 500 GeV, m; = 300 GeV, Af = 300 GeV, // = 200 GeV 
and M2 = 150 GeV with the GUT-hke gaugino mass relations (4). In this case, the 
LSP is bino-like. The spin-independent cross section for the elastic scattering of the 
LSP off proton (aJJ,), the spin-dependent one (cB?), the relic density {fl^h"^) and 
the absolute value of the electron EDM {\de\) are plotted as functions of the phase 
9 Ml in windows (a), (b), (c) and (d), respectively. The sohd line, the dashed hne 
and the dotted line in each window represent the results for tan/? = 10, 30 and 50, 
respectively. In window (c), the region between the two dash-dotted lines corresponds 
to the 2a range of the WMAP constraint 0.0946 < fl^h"^ < 0.1286. In window (d), 
the region below the dash-dotted hne is allowed by the experimental bound for the 
electron EDM \dJ < 1.6 x lO-^'^e • cm. 
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Figure 2: The same as Fig. 1 for tua = 500 GeV, rrif = 800 GeV, Af ^ 800 GeV, // 
= 200 GeV and M2 = 300 GeV with the GUT-hke gaugino mass relations (4). This 
corresponds to mixed LSP with bino-higgsino mixing. 
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Figure 3: The same as Fig. 1 for tua = 500 GeV, = 1 TeV, Af = 500 GeV, /i 
= 700 GeV and M2 = 2.5 TeV with the GUT-hke gaugino mass relations (4). The 
LSP in this case is higgsino-hke. Note that in our analysis on Q^/i^, we did not 
include coannihilation effects. For the higgsino-like LSP, inclusion of this effects will 
significantly reduce fl^h^ so that the lines in window (c) will be greatly modified. 
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rrij = 500 GeV, mj = 2 TeV, A, = 500 GeV, fj. = 2 TeV. = 1.8 TeV, AMSB-like 
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Figure 4: The same as Fig. 1 for uia = 500 GeV, = 2 TeV, Af = 500 GeV, = 
2 TeV and M2 = 1.8 TeV with the AMSB-hke gaugino mass relations (5). The LSP 
in this case is wino-like. 
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nij = 500 GeV, my = 300 GeV, A, = 300 GeV, \fj,\ = 300 GeV, = 150 GeV, GUT- 
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like 




Figure 5: The 9^ dependence of a^^, cr^, ^l^h'^ and |cie| for the same choice of 
rameters as in Fig. 1 (the bino-hke LSP case). 
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Figure 6: The 9^ dependence of a^^, a^, VL^h? and |cie| for the same choice of 
rameters as in Fig. 2 (the mixed LSP case). 



23 



nij = 500 GeV, my = 1 TeV, A, = 500 GeV, = 700 GeV, Mg = 2.5 TeV, GUT-like 
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Figure 7: The Q^^ dependence of c^p, 0"^^ , fi^^/i^ and |cie| for the same choice of 
rameters as in Fig. 3 (the higgsino-hke LSP case). 
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Figure 8: The dependence of c^p, cr^p , VLy}{^ and |cie| for the same choice of pa- 
rameters as in Fig. 4 (the wino-hke LSP case). 



25 



